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Repeating fast radio burst 20201124A originates 
from a magnetar/Be star binary 


F. Y. Wang@ '2™" G. Q. Zhang", Z. G. Dai’? & K. S. Cheng* 


Fast radio bursts (FRBs) are cosmic sources emitting millisecond-duration radio bursts. 
Although several hundreds FRBs have been discovered, their physical nature and central 
engine remain unclear. The variations of Faraday rotation measure and dispersion measure, 
due to local environment, are crucial clues to understanding their physical nature. The recent 
observations on the rotation measure of FRB 20201124A show a significant variation on a day 
time scale. Intriguingly, the oscillation of rotation measure supports that the local contribu- 
tion can change sign, which indicates the magnetic field reversal along the line of sight. Here 
we present a physical model that explains observed characteristics of FRB 20201124A and 
proposes that repeating signal comes from a binary system containing a magnetar and a Be 
star with a decretion disk. When the magnetar approaches the periastron, the propagation of 
radio waves through the disk of the Be star naturally leads to the observed varying rotation 
measure, depolarization, large scattering timescale, and Faraday conversion. This study will 
prompt to search for FRB signals from Be/X-ray binaries. 
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integral of the line-of-sight component of magnetic field 

strength and electron density, and dispersion measure 
(DM)!, defined as the line-of-sight integral over the free electron 
density, of repeating FRBs contain essential information about 
their physical nature and central engine, which are as yet 
unknown!~*, Recently, the rotation measure of FRB 20201124A 
shows dramatically oscillations®. FRB 20201124A was detected 
by the Canadian Hydrogen Intensity Mapping Experiment 
(CHIME)® on 24 November 20207. It has a dispersion measure 
of 413.52 pccm™, larger than the Galactic DM contribution. 
This FRB is extremely active, with burst rate up to 50 per hour?. 
Thanks to its repetition property, it was localized to a galaxy 
named SDSS J050803.48 + 260338.0 at z= 0.098 by Australian 
Square Kilometre Array Pathfinder (ASKAP)®, the Five- 
hundred-meter Aperture Spherical radio Telescope (FAST)°, 
Very Large Array’, and European Very Long Baseline Inter- 
ferometry Network!®. 

FAST detected 1863 independent bursts in a total of 88 hours 
from 1 April to 11 June 2021, covering 1.0 GHz to 1.5 GHz°. 
FAST observations discovered three previously unseen char- 
acteristics of repeating FRBs. First, this FRB shows dramatic 
Faraday rotation measure variations on a day time scale. The 
maximum variation of RM is about 500 rad m~?. A small var- 
iation of RM (about 1.8%) was also found from 20 bursts 
observed by the Effelsberg Radio Telescope!!. Second, it is the 
first repeating FRB shows circular polarization. Third, it is the 
repeating FRB that has the widest mean burst width!?. These 
characteristics cannot be explained by current theoretical models. 
FRB 20201124A can be analogy to the pulsed emissions from the 
Galactic binary system PSR B1259-63/LS 28831, containing PSR 
B1259-63 and the Be star LS 2883. The pulsed emission of PSR 
B1259-63 also shows variations in its RM, DM, scattering time- 
scale and polarized intensity on short time-scales when it 
approaches periastron. It has been argued that a highly magne- 
tized neutron star in an interacting high-mass X-ray binary sys- 
tem can account for all of the observational phenomena of FRB 
20180916B!4. A Be/X-ray binary system, which is composed of a 
neutron star and a Be star with a decretion disk can well explain 
the frequency-dependent active window of FRB 20180916B)». 

In this work, we propose that FRB 20201124A is produced by a 
magnetar orbiting around a Be star companion with a decretion 
disk. The interaction between radio bursts and the disk of the Be 
star can naturally explain the varying RM, depolarization, large 
scattering timescale, and Faraday conversion of FRB 20201124A. 
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Results 

PSR B1259-63/LS 2883. PSR B1259-63/LS 2883 is a binary 
system hosting a radio pulsar and a high-mass Be star LS2883. 
The spin period of the pulsar is 47.76 ms and the spin-down 
luminosity is 8 x 103° erg/s!®. Its distance is 2.6kpc!’. The 
orbital period is 1237 days with eccentricity e = 0.9!°-!8, LS2883 
is an early Be type star that possesses a powerful stellar wind and 
an equatorial decretion disk!?. There is some evidences that the 
Be-star disk is highly inclined to the orbital plane both from 
the pulsed emission?°, the unpulsed emission*!, and timing 
measurements~’. 

This system shows unpulsed emissions in multiple bands, such 
as in radio!®, X-rays*?, GeV24*5, and TeV y-rays?®. The light 
curves show double-peak profiles, which can be attributed to the 
interaction between the pulsar radiation and the disk’. 

Four pulsed radio observations near periastron have been made 
in the year of 19941, 19979, 200078 and 2004?°. When the pulsar 
approaches the disk, the pulsed emissions of PSR B1259-63 show 
large variations in flux, degree of linear polarization, DM and 
RM!3829_ These variations can be explained as the pulsed 
emissions interacting with the disk!33°, The properties of pulsed 
emissions, such as degree of polarization, DM and RM, are 
similar for the four periastron passages”?. 


Rotation measure variation. We attempt to explain the observed 
properties of FRB 20201124A using a physical model for the 
decretion disk of a Be star (see methods, subsection The disk 
model). Our model is shown in Fig. 1. The Be star locates at the 
disk center. For the coordinate system, we set the Be star as the 
origin and fix the x-axis and y-axis on the disk. The orbit has an 
inclination ¢ with respect to the disk. Generally, a magnetar is 
formed by a supernova explosion. A supernova explosion is 
usually asymmetric and gives a large kick to the newly formed 
magnetar*!. Due to the natal kick of the magnetar, the orbital 
plane and the disk plane are misaligned*”. From observations, the 
inclination angle has a wide range in Be star/X-ray binary sys- 
tems, from 25° to 70°°3. In our model, after the companion dies as 
a supernova, a magnetar is formed, leading to FRBs generated in 
the magnetar magnetosphere>*4, which is supported by the FAST 
observations. As the orbital motion of the binary, radio bursts 
from the magnetar pass through different components of the disk 
to reach the observer, which causes the variation of RM. When 
the pathway of radio bursts does not intersect the disk, however, a 
non-evolving RM is expected. 
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Fig. 1 A schematic diagram of the magnetar/Be star binary model. a The Be star locates at the center of the disk. The magnetar is shown as the red point. 
The stellar disk plane (shown as the black dash line) is inclined to the orbital plane (shown as purple dash lines) with the angle ¢ = 2/12. The purple dash 
line is the major axis of the magnetar's orbit. When the pathways of radio bursts pass through the disk, the interaction between bursts and disk can 
reproduce the observed variable rotation measure, depolarization, large scattering timescale, and Faraday conversion. b The face-on view of the Be star's 
disk. The magnetic field shown as arrows is assumed as azimuthal (or toroidal) in the disk. This model predicts that the RM contribution from the disk 
changes sign when the magnetar passes in front of the Be star. The RM variation discovered by FAST supports this scenario. 
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Fig. 2 The fit of RM evolution of FRB 20201124A using the binary model. 
The blue scatter is the mean of observed ARM in day, and the red line is the 
predicted evolution of our model. The blue shadow region is the range of 
observed RM. In the calculation, the disk density is pp = 3 x 10—'4 g cm-3 
close to the star and has a steep-decay index = 4. The orbital period is 
taken to be 80 days with orbital eccentricity e= 0.75. Our model can 
reproduce the main structure of RM evolution. However, there are some 
minor structures can not be reproduced. These minor structures may be 
caused by the clumps in the disk. Source data are provided as a Source 
Data file. 
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For the binary system PSR B1259-63/LS 2883, significant RM 
changes were observed both in magnitude and in sign before and 
after the 2004 periastron passage*?. Compared to the Galactic 
PSR B1259-63, the RM variation of FRB 20201124A is small. A 
thin disk is considered. We assume the density of disk at the 
stellar surface pp = 3 x10~!4gcm~3 and B=4, which are the 
typical value for disks of Be stars*>. The orbital period is assumed 
to be 80 days with orbital eccentricity e= 0.75. Using the disk 
model in equation (1) and the magnetic field profile in equation 
(2), we can reproduce the overall RM variation with time. For the 
spherical coordinate origin at the Be star, the observer’s direction 
is @=3.08 rad and ¢=2.18 rad. The derived RM variation is 
shown as the red line with magnetic field By = 10 G at the stellar 
surface in Fig. 2. This magnetic field strength is consistent with 
those derived from RM variation and Faraday conversion. The 
blue points is the mean observed RM variation ARM = RM-RM, 
in a day. RM, is taken to be constant, as observed after MJD 
59350. The blue region shows the range of ARM in each day. 

Some dips and peaks of ARM are not well reproduced. The 
possible reason is that the disk is very clumpy both in density 
and/or magnetic field. From observations, a clumpy disk was 
required to explain the multi-wavelength observations of the 
Galactic binary system PSR B1259-63/LS 2883°°?7. The clumpy 
disk is also supported by numerical simulations. For example, a 
two-armed spiral density enhancement caused by the tidal 
interaction was found using three-dimensional smoothed particle 
hydrodynamics simulation*®. From observations, a stellar wind 
from a massive star is clumpy, with the filling factor of f, = 0.1°°. 
At present, there are no constraints on disk clumps. For 
simplicity, we assume that they are similar to the clumps in a 
stellar wind, without magnetic field clumps. The density inside 
the clumps is 1/f, larger that of the smooth disk. Because the 
interspace between clumps is not empty, the clumping factor can 
be expressed as f, = (n2)/(n,)?. The spatial distribution of 
clumps remains unknown. Using ad hoc clumps to fit the RM 
variation is less productive. So we discuss the effect of clumps on 
RM quantitatively. For the RM peak near MJD 59340, the 
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Fig. 3 The DM and free-free absorption optical depth z contributed by the 
disk. The same parameters as Fig. 2 are used. The DM variation shown as 
the blue line is small, which is consistent with the FAST observation. The 
optical depth z due to free-free absorption by the disk denoted as the red 
line is much less than 1. So the existence of the disk has no effect on the 
detection of radio bursts. 


observed value is about ten times larger than the predication of 
the smooth disk. This implies the value of f, is about 0.1 if only 
density clumps are considered, being similar to those in a 
stellar wind. 

The observations show that some bursts have low polarization 
before MJD 59350. The fluctuations of the electron density in the 
disk can cause differential Faraday rotation along different paths, 
which can account for the observed burst depolarization (see 
methods, subsection Depolarization in the disk). The mean burst 
width of FRB 20201124A is larger than the widths of other known 
repeating FRBs. The fluctuations of the electron density also can 
cause scatter broadening of the pulses (see methods, subsection 
Scattering timescale). 

The disk also affects the DM and detection of radio bursts. In 
Fig. 3, we show the DM and the optical depth due to the free-free 
absorption contributed by the disk. The DM change shown as the 
blue line is small. The maximal change is 0.12 pccm~>. FAST 
found no significant DM variation in a 95%-confidence-level®, 
with an upper limit ADM <4.9 pccm~?. Therefore, the model 
prediction is consistent with FAST observations. Next, we 
consider the free-free absorption as a function of orbital phase. 
The optical depth 7 due to free-free absorption shown as the red 
line is less than unity. The DM and free-free absorption caused by 
the stellar wind can be safely ignored (see methods, subsection 
The effect of a stellar wind). The presence of disk clumps also 
affects the interaction with the radio bursts. The DM and free-free 
absorption are enhanced by a factor of 1/f,= 10 with respect to 
the case of a smooth disk*°. From Fig. 3, the maximum DM 
contribution is 4pccm~? in the clumps, which is less than the 
observed DM variation®. After considering the disk clump, the 
optical depth to free-free absorption is less than 1. So the clumpy 
disk does not affect the detection of radio bursts. For FRBs, the 
induced Compton scattering should be considered*!. The induced 
Compton scattering due to the disk is also estimated to be small, 
and does not affect FRB detection. Therefore, the disk is 
transparent to radio emission, and has no effect on the detection 
of FRBs. 


The circular polarization from Faraday conversion. For non- 
repeating FRBs, variations in both linear and circular polarization 
have been observed**. However, previous observations show that 
repeating FRBs only have high degrees of linear polarization. FRB 
20201124<A is the first one showing high circular polarization. For 


NATURE COMMUNICATIONS | (2022)13:4382 | https://doi.org/10.1038/s41467-022-31923-y | www.nature.com/naturecommunications 3 


ARTICLE 


NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-022-31923-y 


0.100 F 3 


0.010 F 3 


0.001 F 3 


Faraday conversion coefficient pg 


1000 


minimum Lorentz factor Yin 


Fig. 4 Faraday conversion coefficient pq distribution as a function of 
minimum Lorentz factor y,,,,- An electron column density of N.=1pc cm—3 
is assumed. The plot uses vy = 1.25 GHz, 0= 7/4 and p= 2.3. Faraday 
conversion only takes effect for y,,;, < 200. 


the Galactic binary system PSR B1259-63/LS 2883, the linear and 
circular polarizations varied rapidly in an hour timescale near the 
periastron2*. Physically, besides Faraday rotation, the propagation 
of FRBs through a magneto-ionic environment can also lead to 
Faraday conversion. In addition, a fraction of circular polariza- 
tion of FRB 20201124A is produced by radiation mechanism». 
There are two mechanisms to generate Faraday conversion. The 
first case is that radio signals propagate into a relativistic 
plasma*?. The second case is that the magnetic field quasi- 
perpendicular to the wave vector is large, which had been used to 
interpret the circular polarization observed in solar radio bursts*+. 
For the first case, termination shocks can accelerate electron and 
positron pairs from the magnetar wind to relativistic speeds. 
Because electrons and positrons contribute with the opposite sign 
to the circularly polarized component, the different distributions 
of electrons and positrons are required to generate circular 
polarization. In our model, a weak magnetar wind is considered. 
Therefore, the termination shock may be absent. For the binary 
system PSR B1259-63/LS 2883, the minimum Lorentz factor of 
the shocked electrons is estimated as*° 10°. In Fig. 4, the Faraday 
conversion coefficient as a function of minimum Lorentz factor 
Ymin 18 Shown. From this figure, we can see that the Faraday 
conversion coefficient is found to approach zero when y,,,, is 
about 200. From the accurate calculations, the exact Faraday 
conversion is found to approach zero with the increasing electron 
energy“, because ultra-relativistic electrons hardly interact with a 
radiation field. When the magnetar is in front of the Be star, FRB 
signals cannot interact with the plasma shocked by the termi- 
nation shock, leading to no conversion. So, the Faraday conver- 
sion due to relativistic plasma is hard to account for the observed 
circular polarization (see methods, subsection Origin of the cir- 
cular polarization). 

Next we use the magnetic field quasi-perpendicular to the wave 
vector to explain the origin of circular polarization. When the 
cyclotron frequency is approaching the observation frequency, 
the Faraday conversion will take effect. From the FAST 
observations, the required magnetic field B, is about 7 Gauss 
for Faraday conversion®, which is the same order of the parallel 
magnetic field By. This mechanism naturally explains the circular 
polarization of FRB 20201124A. FAST observations indicate that 
the Faraday conversion occurs in the RM variation stage (e.g., 
bursts 779 and 926). When the RM variation stops, the circularly 
polarized bursts lack qusai-periodic structures (e.g. Burst 1472). 
There is no Faraday conversion. This is consistent with the 
prediction of our model. For the above orbital parameters, the 


observed Faraday conversion appears when bursts’ paths 
penetrate the disk. When the pathway of radio bursts does not 
intersect the disk, the RM variation stops. Meanwhile, there is no 
perpendicular magnetic field required for Faraday conversion. 


Candidates of a magnetar/Be star binary. The total energy 
emitted by radio bursts during FAST observation is about 
3 x 104! erg in 88h°. For dipolar magnetic field Bg = 10) G, the 
total magnetic energy of a neutron star (NS) is E;ag= 3 x 108 
erg. In order to explain the energy budget of FRB 20201124A, 
about 3% NS magnetic energy is transferred to radio bursts in 
88h, which is impossible. Moreover, the Galactic FRB 
20200428 suggests FRBs are accompanied by X-ray bursts 
(XRBs). The energy ratio between a radio burst and its associated 
XRB is about 10~> 4°. If this value is valid for FRB 20201124A, 
the energy of XRBs is Ex =3 x 104!/10~> ~ 104° erg. For Emag > 
Ex, the magnetic field B, is larger than 101° G. So it should be a 
magnetar. The bursting activity may be related to the crustal 
magnetic energy of young magnetars*/. If the converted magnetic 
energy comes from the toroidal component, it is at least two 
orders of magnitude larger than the dipolar magnetic field of the 
neutron star. From numerical simulations*® and observations*’, 
the toroidal component could be ten times larger than the dipolar 
component. So, if the required magnetic field is toroidal 
(>10!° G), the dipolar magnetic field is larger than 10!4G, which 
also supports a magnetar. Therefore, if the central engine is a 
compact stellar object, it should be a magentar from energy budget. 

There are two candidates of a magnetar/Be star binary. The 
first one is LS I + 61° 303°°. The second questionably one is SGR 
0755-2933, which is also referred to as 2SXPS J075542.5- 
293353°!. LS 1+ 61° 303 is a Be gamma-ray binary, with period 
about 26.496 d and a high eccentric orbit (e = 0.537)°*. The mass 
of the Be star in this binary is about 13 Mo. It has been found that 
some magnetar-like bursts are possibly from this binary°®?. 
More recently, transient radio pulsed emissions from its direction 
were detected by FAST°4, with a period of 269 ms. This supports 
the presence of a rapidly rotating neutron star in this binary>4. So 
the magnetar/Be binary used in this paper is close to LS I+ 61° 
303. They have similar orbital period and eccentricity. Searching 
for FRBs from LS I+ 61° 303 is attractive. From the FAST 
observation, it has been found a rapid change of the plasma 
properties along the line of sight, which may be caused by the 
variation of the stellar wind properties*+. Because the brightness 
temperatures of FRBs are much higher than that of pulsations 
of neutron star°°, the induced Compton scattering may be 
significant4! >. 

However, no FRB-like signal from the two systems has been 
reported. The possible reasons are as follows. First, FRBs are rare 
phenomena. They may be highly beamed>”*®, and have narrow 
spectra, making them hard to discover. For example, more than 
20 magnetars have been discovered in our Galaxy. Only FRB 
200428 emitted from SGR 1935+2154 has been discovered??. A 
bright X-ray burst is associated with FRB 200428, which has a 
strikingly different spectrum from other X-ray bursts®!?, FAST 
observed SGR 1935+2154 for eight hours in its active phase°”. 29 
X-ray bursts are detected, but no FRBs coincident with the bursts. 
Second, the high brightness temperatures of FRBs require that the 
emission process mechanism must be coherent!34, The physical 
conditions required to produce coherent radiation in magnetars 
may be hard to satisfy’. Third, the burst rate of FRBs deviates 
from Poissonian and varies significantly®*-®. For FRB 
20201124A, FAST found a sudden quenching of burst activity°. 
If similar active properties of FRBs are produced in magnetar/Be 
star binaries, enough and proper-time observations should be 
performed to search FRB signals from them. 
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Fig. 5 The fit of RM evolution for FRB 20190520B using the binary 
model. The blue points are source RM calculated from RM,, x (1+ 25)”, 
where RMop is the observed value and Zop = 0.241 is the redshift. The red 
line is the fit from the binary model. The disk density is pp = 8 x 10-8 
g/cm with a steep-decay index of — 3. The mass of the Be star is assumed 
to be 30Mg with radius R, =10Ro. The magnetic field at the stellar surface 
is 100 G. The orbital period is taken as 600 d with orbital eccentricity 
e= 0.75. Source data are provided as a Source Data file. 


Discussion 

Considering the interaction between a Be star disk and FRBs, we 
can naturally explain the mysterious features of FRB 20201124A. 
This model can well produce the variation of RM with a magnetic 
field of 10 Gauss at the stellar surface. This magnetic field and the 
strong electron density fluctuations can generate the observed 
depolarization of the FRBs by differential Faraday rotation along 
different paths. The magnetic field perpendicular to the wave 
vector in the disk, having the same order of the magnetic field 
along the line of sight, can account for the observed circular 
polarization. 

Recent observations show that the RM of FRB 20190520B is 
rapidly varying with two sign changes®®°’, with similar magni- 
tude as that of PSR B1259-63/LS 2883. The sinusoidal-like form 
of RM variation emerges. A dense magnetized screen near the 
FRB source can account for the RM _ variation and 
depolarization®. In our model, the distance of the screen can be 
of order the separation of the two stars. For orbital period 
P=600 d, our model can well explain the RM variation of FRB 
20190520B, which is shown in Fig. 5. This figure gives the total 
RM in the source frame. The RM contributed besides the disk is 
assumed to be RM,(1 + Z,,) = —10301.92 radm~?, where 
RM,» = — 6700 rad m~? is the first observed RM by Green Bank 
Telescope (GBT) at MJD 59106°° and z,, = 0.241 is the redshift 


of FRB 20190520B. The semimajor axis can be derived as a= 


(GMP? /2n)'? = 8.5 AU, where Mr = 30Mo is the total mass of 
the binary. Our model predicts the RM evolution would be quasi- 
periodic. If a large number of RM detections spanning a long 
timescale are accumulated, the orbital period could be derived 
from these RM data. 


Methods 
The disk model. There are some evidences for a disk around a Be star? and 
especially for LS 2883!°. The disk opening angle is found to be less than 20°. The 
disk size ranges from a few stellar radius to a few tens stellar radius from inter- 
ferometry observations**. The density of disk can be well described by a radial 
power-law form, with surface density between about 10~}3 to a few times 
10—!0 g cm~3 35. From the DM and RM evolution of PSR B1259-63/LS 2883, 
Po = 10~!6 gcm~3 was used?9, 

The classical formation channel of a Be/X-ray binary (or a magnetar/Be star 
binary) is as follows. In an interacting binary, the B type star accretes matter from 


the companion star. Therefore, it is spun up to nearly critical velocity, leading to 
the formation of a decretion disk, which is supported by observations”. At this 
stage, the orbit and disk are almost in the same plane”!. In a magnetar/Be star 
binary, the companion is a magnetar, which is usually formed by a supernova 
explosion. The process of supernova explosions is asymmetric and gives a very 
large kick to the newly formed magnetar*!-”*. Due to the kick of the magnetar, the 
orbital plane and the disk is misaligned. Besides, the kick also makes the orbit be 
eccentric>*79, 

From observations, some Be/X-ray binaries have been found to be high 
misaligned>>. For the Galactic binary PSR B1259-63/LS 2883, the inclination angle 
between the orbital plane and the disk plane are measured using several methods 
from observations. The inclination angle is found to be 10° to 40° by fitting to the 
observed variations in DM and RM of PSR B1259-63°. Using the X-ray data of 
PSR B1259-63, the derived inclination angel is about 70°74. The inclination angle of 
PSR B1259-63/LS 2883 system is about 25° from the optical spectra of LS 288375. 
Using the 23 yr of pulsar timing data, the inclination angel is found to be 35°!8. 

We consider the rotationally supported, geometrically thin disk model with 
isothermal hydrogen gas. The density profile of disk can be described as*° 


-B 
pr, Z) = Po (=) exp[—z*/H*(r)], Q) 


where r is the distance from the star, pp is the disk density at the stellar radius R,, 8 
is the density slope, z is the height, and H(r) is the scale height of the disk. 
Observational properties of Be stars are well described by this disk model and the 
density slope B is from 2 to 495. 

The magnetic field, B, in the disk is poorly known. We assume that it is 
azimuthal and has a form 


B(r) = B,(R,/7), (2) 


with Bo is the magnetic field strength at the stellar surface. This model can explain 
the RM variation of PSR B1259-63/LS 2883!°. Given the above model, we can 
calculate the RM as a function of orbital phase. The RM contributed by the free 
electrons in the disk is 


RM = 8.1x 10° [ n,B,dlradm~, (3) 
LOS 


where n, is the free electron density in unit of cm~3, B) is the magnetic field along 
the line of sight in unit of Gauss, andlis in unit of parsec. For the azimuthal 
magnetic field described above, we expect the RM contributed by the disk to 
change sign when the magnetar passes in front of the Be star, which is supported by 
FAST data. 


The effect of tidal forces on the disk. The tidal force of the magnetar could affect 
the disk’®. The process of disk truncation in Be/X-ray binaries was first discussed in 
1997 from the orbital period-Ha@ equivalent width diagram’’. From observations, 
26 Be stars show disk truncation signature’®. Among them, six Be stars are known 
to have close binary companions in circular orbits’®. In the classical viscous disk 
model”, the criterion for disk truncation can be written as 


dg als Tres <0, (4) 


where T,;, and T,,, are the viscous torque and the resonant torque, respectively. 


Using the standard formulae of Tyi, and T',.,°°, Ty, + Tes $0 is met for a small 
value of « (viscosity parameter). The drift time-scale is Taig = Ar/v,, where v, is the 
radial velocity, and Ar is the gap size between the truncation radius and the radius 
where the gravity by the magnetar begins to dominate. The spread of the disk can 
make the truncation inefficient for a system with Tgrig/Porb = Ar/0.1¢,Porp < 1. For a 
typical value a = 0.1, Tari¢/Porb = 0-1 in our model, which supports the truncation 
effect is inefficient. This conclusion was also confirmed by numerical simulations. 
For misaligned Be/X-ray binary systems, smoothed particle hydrodynamics 
simulations show that Be star disks are not truncated in highly eccentric binaries>®. 
In addition, the disk may be warped by the tidal force*®, leading to non-smooth 
distributions of density and magnetic field in the disk. This is the possible reason to 
cause the clumps in the disk. The spread of disk could cause the smaller density of 
the disk’s outer part, compared to last periastron. Therefore, the variation of RM 
will become smaller in magnitude over time. In our model, the Be/magnetar binary 
has a highly misaligned disk, eccentric orbit and ignorable disk truncation. These 
three requirements are self-consistent. 

Besides the tidal force, the magnetar wind may also affect the disk. In our 
model, the radio bursts can be powered by the magnetic-energy release of the 
magnetar. The magnetar wind may be weak, similar to that of Galactic magnetars 
(ie, SGR 1935+2154). Strong magnetar winds can produce a bright, compact 
persistent radio source associated with FRBs*!,8?, which has been observed in FRB 
20121102A3 and FRB 20190520B*4. However, there is no compact persistent radio 
emission associated FRB 20201124A!°. The extended radio emission around FRB 
20201124A is attributed to star-formation activity®?®, 


Fitting RM using the binary model. Compared to the Galactic PSR B1259-63, the 
RM variation of FRB 20201124A is small. We consider a thin disk by assuming 
Po=3x10-!4gcm-3, and B =4. The radius and mass of Be star is taken to be 
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Ry, =5Ro and M, = 8Mg. Typical velocities of the plasma in the disk are®® vq = 
150—300 km s~!. The vertical scale height is derived from 


7 \'P 
= 5 
Hn=c (=) r, (5) 


where c, is the sound speed. For an isothermal gas, the sound speed is 


c, = (kT /um a! > where T is the isothermal temperature, y is the mean molecular 
weight, and my, is the mass of hydrogen. We assume ps = 1 and T= 10‘ K. In order 
to prevent that the disk is significantly truncated by tidal forces, we consider a large 
eccentric orbit with e = 0.75. In our model, the RM variation is mainly affected by 
the orbital motion. We assume the orbit period is P= 80 days. Then the semi- 
major axis can be derived from 


G(M, + Mys)P?]"? 

ne aoe (6) 
4n? 
where Mys = 1.4MoQ is the magnetar mass. To match the observed RM evolution, 
we assume that the inclination angle is ¢ = 1/12 and the observer's direction is 
0 = 3.08 rad and ¢ = 2.18 rad. We assume the magnetic field Bp = 10 G. From 
observations, the surface magnetic field of a Be star has been measured. The 
magnetic field in the longitudinal direction is between 2 and 300 Gauss for 15 Be 
stars®”. From a large sample of 60 Be stars, the maximal surface magnetic field 4000 
Gauss was derived*®, Therefore, the value of By = 10 G is well in the allowed range. 
The DM contributed by the disk is 


DM = j n,dl, (7) 
LOS 


where n, is the free electron density along the line of sight. In the observations, 
FAST can detect FRB signals. Therefore, the optical depth 


d 
T= 8.2x10?T Sy?! | nedl (8) 
0 
of free-free absorption must be less than unity. In the above equation, T is the 
temperature in Kelvin, v is the observing frequency in GHz, and d is the distance to 
the magnetar (in parsec). The same parameters for the disk is used to calculate the 
DM and optical depth. 

The induced Compton scattering is crucial for high brightness temperature 
sources. In our model, the scattering occurs in the nearest vicinity of the source. 
For a single narrow beam of uniform brightness and solid angle AQ, the effective 
optical depth is#!>° 


kT 
2 (AO )ty, (9) 


Ting = 
in m,c 


where T; is the brightness temperature, and tp=1,(r)o;r is the Thompson 
scattering optical depth. For radio bursts with duration w and flux S occurring at 
redshift z= 0.098, the brightness temperature is 


_ o72 Bi sy 30 s BN? WN d ° 0 
Ty, = Sd2/2k(uwy ~ 6.6x 10 «(F) (Gs) (—) (same) - (10) 


The radiation subtends a solid angle AQ = (cw)?/r3, = 10~!° with the distance 
rg=50Ro between the magnetar and the disk. If we assume radio bursts travel at 
distance 2R, away from the Be star, using the base density 5 x 10—!4 gcm~? of the 
disk and equation (9), we can estimate T;yq is 0.025, which is much less than one. 
Therefore, radio bursts can escape from the disk. 


The effect of a stellar wind. A massive star usually has a large wind mass loss, 
which strongly depends on the stellar effective temperature®’. The effect of a stellar 
wind on DM and free-free absorption should be considered. For an isotropic wind, 
the number density of the wind is given by 


Pont = frond BO: (11) 


where r is the radial distance from the star. The stellar surface density is py \;.g = 
M/4nR2v,, 
observations, the mass loss rate of Be stars has been estimated to be 
10-!!10-8Mo yr7!. A typical value powind = 4-18 x 10—!7 gcm~? is adopted, 
which corresponds to the mass loss rate 3 x 10~-!9Mg yr—!. The ratio of the mass 
flux density in the disk over the stellar wind is 


Poa 
Po,wind Vw 
which is well in range of 30-100%. 
Because of adiabatic cooling, the wind temperature decreases slowly with radial 


distance from the star. The temperature of the stellar wind follows a power-law 
evolution due to adiabatic cooling? 


T wind — To wina(t/R,) 2 (13) 
In our calculation, we use the effective temperature of star Towing = 3 X 104K 
and $, = 2/3. In this case, the maximum DM,, and 1,, contributed by the stellar 
wind are 0.54 pc cm~3 and 0.005, respectively. Based on these results, the 


with mass loss rate M and wind velocity v,,=3 x 108cms~!. From 
90,91 


~ 80, (12) 


contribution of the stellar wind can be safely ignored. We are not aware of 
magnetic field measurements in the wind of a massive star. For the Galactic PSR 
B1259-63, the RM variation cannot be explained by the contribution from the 
stellar wind*°. The geometry of the magnetic field in the stellar wind may be 
toroidal>°.94, Therefore, the radio burst propagates orthogonally to the toroidal 
magnetic field in a significant fraction of the orbit, which contradicts the observed 
RM evolution. 

Below, we compare the magnetic field in the pulsar wind and in stellar disk. 
Since the magnetic field is dipolar inside the magnetar light cylinder, 
Tuc =5 x 10°P,,, cm, and toroidal (« r~!) in the wind, the magnetic field strength 
at a distance r is 


_;/5x10? cm 
Byina = 230B,4P, a r ) G, (14) 


for the period of magnetar P,,, = 2 s. We make the comparison at the same radius 
r=2R, = 10Ro. Bywind = 0.9 G is less than the magnetic field in the disk. Therefore, 
the magnetic field in the disk dominates. 


Depolarization in the disk. The observed pulse depolarization is attributed to 
random fluctuations in the magnetic field and/or electron density in the disk, 
which can cause differential Faraday rotation along different paths of the light 
rays”>. If the Faraday dispersion function is well described by a Gaussian dis- 
tribution with variance o” = k?(5(RM)*), the measured polarized intensity can be 
parameterized as 


P()’) = P; exp[—2k? (5(RM)*)A“], (15) 


where P; is the intrinsic polarized intensity, A is the observing wavelength, k = 0.81 

; 1/2, : : eee 

is a constant, and (6(RM)’) ? is the variance in ARM. For simplicity, we assume 
1/2, eee : 

that (6(RM)") ? arises from fluctuations in free electron density n,. Furthermore, 


we assume that the variance of electron density n, satisfies (6n,2)1/ 2 & n, in the 
disk>°, which is similar to the small-scale turbulence in the interstellar medium. 


Letting (o(RM)?) = fARM, the stochastic Faraday rotation leads to the depo- 
larization when |ARM| > f—!A~2. From the FAST observation, the value of |ARM| is 
about 200 rad m~? at frequency 1.25 GHz. Therefore, f is about 0.1. We conclude 
that stochastic Faraday rotation may be responsible for the observed 
depolarization. 


Scattering timescale. From the CHIME observation, the largest scattering time- 
scale of FRB 20201124A is 14 ms’. The fluctuations of the electron density in the 
disk can cause scatter broadening of the pulses. The electron plasma with density 
fluctuations of dn, and spatial scale a produces a typical scattering angle at 
wavelength A of 


1 /L\2 
0, = 5 () rd’ 6n,, (16) 
where fo is the classical electron radius and thickness of the screen L is assumed to 
be roughly equal to its distance d from the magnetar”®. If the scattering is due to a 
thin screen at a distance r, from the magnetar, the scattering angle 0, is (2T,c/ ri! A 
with the scattering time t,. We use the burst 20210405B with scattering time 

T, = 12.17 ms detected by CHIME’ to estimate the spatial scale a. The distance 
between the burst and the screen is about 42 solar radius. Therefore, the scattering 
angle @, is about 0.9°. From equation (16), the density fluctuation of electron 
plasma in the disk is dn, = 80a!/? cm}. If dn, is set to be equal to the total electron 
density 1, in the scattering region, from the disk model, then dn, = 3.9 x 10°cm~3 
is found, implying that a = 2.4 x 109 cm. 


Origin of circular polarization. FRB 20201124<A is the first repeating FRB with 
significant circular polarization as discovered by FAST? and Parkes”’. Two possible 
mechanisms can cause Faraday conversion: linearly polarized light is converted to 
circularly polarized one and vice versa. The polarization measurement of FRB 
121102 has been used to constrain the local magneto-ionic environment using 
Faraday conversion?®9, 

First, we consider that the circular polarization arises from the propagation 
through mildly relativistic plasma. For the PSR B1259-63/LS 2883, X-ray, GeV 
gamma-ray and Tev gamma-ray radiations have been observed. The natural 
physical mechanism for these high-energy emissions are as follows. An interaction 
between the relativistic wind of the pulsar and the low-velocity wind from the Be 
star can form a termination shock. The electron and positron pairs from the pulsar 
wind are accelerated to relativistic velocities by the terminal shock and emit 
broadband non-thermal emission through synchrotron radiation and inverse 
Compton scattering. In our model, a termination shock may not exist, due to the 
weak magnetar wind. 

If only the relativistic plasma is considered, the change of polarization angle 
is 

1 3 
Ag => RRM)’, (17) 
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where RRM is the relativistic rotation measure (RRM), and A is the observed 
wavelength. We assume a power-law distribution for accelerated particles 


NY) = Noy?» Ymnin $$ Ymnax> es) 
where y,,;, and y,,,, are the minimum and maximum Lorentz factors, respectively. For 


the power-law index p # 1, the normalization factor is Ny = n,(p — Noy? - ym? me 
with the number density n, of the relativistic particles. The definition of RRM is!° 
4 


c po + |B sin 8) You 


= 19 
4m egmict p — 2 i, 


The fraction of circular polarization is Il, = RRMA? for the relativistic plasma 
dominated case. For relativistic plasma, the definition of DM, is DMx D!°!, with 
Doppler factor D = [y(1 — v/ccos 6] ', where 0 is the angle between the line of sight 
and the direction of plasma motion. The RRM effect should require 

n, <K ny, In(v/f y,)- 

For the 1994 periastron passage of PSR B1259-63, the changes of DM and RM 
cannot be explained in the termination shock model*°. The observed changes of 
DM in four periastron passages of PSR B1259-63 are roughly consistent with the 
contributions from the disk and wind of Be star, with the disk component 
dominated!?, Therefore, the cold plasma dominates the wave properties. The 
RRM effect is very weak. When the magnetar is in front of the Be star, radio bursts 
cannot interact with the plasma shocked by the termination shock. Therefore, there 
is no Faraday conversion, which is disfavored by the FAST observations. 

Numerically, the approximated Faraday conversion coefficient is*3!03 


2 goo? —1 (f,sind\* 2 
Po = 85x 107 = | P= (fesin S, Ch 
p—2)|\f, sin Oyun yoe w w 


where w = 2n1, f, = eB/m,c is the cyclotron frequency, and 0 is the angle between 
wave vector and the line of sight. This approximation approaches zero if y,,,,,, = 200. 
In order to explain the high-energy emission from the Galactic binary system PSR 
B1259-63/LS 2883, an extreme large y,,;, = 10° is required*>!°. Therefore, the 
relativistic plasma hardly generates the observed circular polarization, unless the 
minimum Lorentz factor is less than 100. 

Next we consider the Faraday conversion occurring when the magnetic field is 
perpendicular to the wave propagation direction. The required magnetic field for 
Faraday conversion is about 7 G®, which is the same order of the magnetic field 
strength along the line of sight. The Be star disk can provide the required vertical 
magnetic field. This is the most natural way to explain the production of circular 
polarization. 
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